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ABSTRACT 

Context. Ram pressure stripping of galaxies in clusters can yield gas deficient disks. Previous numerical simulations 
based on various approaches suggested that, except for near edge-on disk orientations, the amount of stripping depends 
very little on the inclination angle. 

Aims. Following our previous numerical and analytical study of face-on stripping, we extend the set of parameters 
with the disk tilt angle and explore in detail the effects of the ram pressure on the interstellar content (ISM) of tilted 
galaxies that orbit in various environments of clusters, with compact or extended distributions of the intra-cluster 
medium (ICM). We further study how results of numerical simulations could be estimated analytically. To isolate the 
effect of inclination, galaxies on strictly radial orbits are considered. 

Methods. A grid of numerical simulations with varying parameters is produced using the tree/SPH code GADGET with 
a modified method for calculating the ISM-ICM interaction. These SPH calculations extend the set of existing results 
obtained from different codes using various numerical techniques. 

Results. The simulations confirm the general trend of less stripping at orientations close to edge-on. The dependence 
on the disk tilt angle is more pronounced for compact ICM distributions, however it almost vanishes for strong ram 
pressure pulses. Although various hydrodynamical effects are present in the ISM-ICM interaction, the main quantitative 
stripping results appear to be roughly consistent with a simple scenario of momentum transfer from the encountered 
ICM. This behavior can also be found in previous simulations. To reproduce the numerical results we propose a fitting 
formula depending on the disk tilt angle and on the column density of the encountered ICM. Such a dependence is 
superior to that on the peak ram pressure used in previous simple estimates. 

Key words. Galaxies: clusters: general - Galaxies: evolution - Galaxies: intergalactic medium - Galaxies: ISM - Galaxies: 
structure 
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■ Ram pressure of the inter-galactic medium in the dense en- 

' vironments of galaxy clusters and groups can significantly 
affect the interstellar content of galaxies. Together with When the galaxy plunges into the intracluster medium 
tidal stripping and gas outflows due to galaxy mergers, (ICM), a complex hydrodynamical interaction of the in- 
ram pressure stripping is an efficient mechanism remov- terstellar matter (ISM) and ICM takes place. Several 
ing their gas component. This quenches the star forma- phases of ISM-ICM int eraction were recognized by 
tion. Local enhancement of the gas density by the ram iRoediger fe Briiggenl (|2006[ ): the initial ram pushing of the 
pressure, however, can increase the star formation rate in ISM out of the disk, followed by stripping of elements 
the central disk regions an d in the tail of stripped mate- accelerated to the escape velocity, in a background sup- 
rial (|Kronberger et al.||2008j). Current ob servations showing ported by continuous stripping. A simple estimate of the 
gaseous tails and extra-planar gas (e.g. |Crowl et al]|2005t stripping amount suggests that if the ram pressure of the 
IChung et aLll2007l). and a number of numerical simul ations wln d exceeds the local gravitational restoring force per 



(e.g. |Vollmer et al.| |2001|; |Roediger fc Bruggen| |2007|) con- un j t ar ea, the gas is releas ed from the potential of the 



firm the importance of gas stripping in the evolution of galaxy (|Gunn fc GottJ I1972L GG72). This approximation 

galaxies. Although some gas tails could be due to tidal wor ks we ll for face-on cases with constant long-lasting flows 

inter actions between galaxies in the cluster (jChung et all f the ICM. However, in real situations, galaxies experience 

I2007D, these are rapidly destroyed by the ram pressure a time- varying ram pressure due to their evolving orbital 

velocity and radially changing ICM density. An approxi- 

Send offprint requests to: P. Jachym mative criterion where the stripping is proportional to the 
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momentum of the ICM encounter ed along the orbit was 
introduced by I Jachym et al.l (|2007f) . 

A long-term ram pressure corresponding to the over- 
all cluster ICM density profile can be accompanied by 
short pressure spikes from local clumps in the ICM dis- 
tribution or wakes cause d by previous stripping events. 
iTonnesen fc Brvanl (|2008l ) studied the role of ICM overden- 
sities and their infalling motions, and suggested that even 
small substructures in the ICM may evoke a stronger strip- 
ping than that corresponding to the given position in the 
cluster. In short ICM-ISM encounters, not all the ISM is ac- 
celerated to t he escape speed and its re-a ccretion can occur. 
Very recently, ITonnesen fc Brvanl ([2009) took into account 
the multi-phase nature of the ISM, and fou nd different re- 
sults, according to the dumpiness of the gas. lVollmer et all 
( 2009f) argue that ram pressure is more efficient on ionized 
gas, and they see evidence of differential ram pressure effi- 
ciency according to the nature of the multiphase gas. 

The question of how the orientation of the galaxy rel- 
ative to the orbital direction influenc es the final strip- 
ping was touched on by many authors (iQuilis et al.|[2 000: 
Vollmer et al.ll200ltlSdiurz fc Struckll200lHHidaka fc Sofuel 
2002; lMarcolini et al.|[2"003tlRoediger fc Bruggenl2006fh For 
tilted galaxies the effect of the hydrodynamical shield- 
ing is more pronounced in comparison to the face-on case 
and yields weaker stripping. Moreover, due to galactic 
rotation, near edge-on ram pressures evoke asymmetri c 
shapes of the ISM disks (e.g. iRoediger fc Briiggenl [20061. 
From previous simulations, only close-to edge-on orien- 
tations were ident ified as less effic i ent in the stripping 
than face-on cases (|Quilis et al.ll2000l : IRoediger fc Bruggenl 

2006) . The tilted stripping was found to be delayed com- 
pared to_Jacfr : on_(^^lIme^^^] [200l|) , and takes a longer 
time ([Farouki fc Shapirolll980t )7Other features, like tails of 
stripped material formin g on the disk side where the gas 
rotates against the wind (jPhookun fc M undv 1991), or an- 
nealing of inclined galaxies by a s tronger loss of angular mo- 
mentum dSchulz fc Struck! l200lh . were noticed. Mimicking 
the hydrodynamical shielding of tilted galaxies w ith a 
simple geometrical prescription, IVollmer et al.l (|200ft ) per- 
formed a set of sticky-parti cle simulations with time- 
varyi ng ram pressure peaks. IRoediger fc Bruggenl (|2006t 

2007) . ' using hydrodynamical adaptive-mesh-refinement 
simulations, introduced models of galaxies experiencing a 
constant ICM flow in a periodic boundary box, or orbiting 
on realistic orbits in clusters. 

What is the role of the galactic disk tilt angle in var- 
ious distributions of the ICM, and could a simple depen- 
dence between the stripping amount and the tilt angle be 
found from simulations? To answer these questions we per- 
form in this paper an extensive grid of Smoothed particle 
hydrodynamics (SPH) simulations with time-varying large 
or small ram pressure peaks of generally tilted ram pres- 
sure stripping events. This completes other existing sim- 
ulations based on different numerical approaches (AMR, 
sticky-particles etc.) a nd follows our simu lations on the 
face-on study made by I Jachym et al.l (|2007h . Rather than 
describing in detail the process of stripping itself (that shall 
be studied in a forthcoming paper), we focus on the quanti- 
tative outcomes of the tilted ISM-ICM interactions. To iso- 
late the effect of the disk inclination, we deal with galaxies 
on strictly radial orbits freely falling from the cluster out- 
skirts. This condition however does not imply a strong lim- 
itation since ram pressure profiles experienced by galaxies 



along slightly elliptical orbits may be modeled with strictly- 
radial ones through an ICM distribution with a lower value 
of the central density and a lower concentration towards 
the center. With the results of the simulation grid in hand, 
we use a simplified test particle model to separate the vari- 
ous effects taking place in the ISM-ICM interactions, and to 
search for a simple dependence that would fit the numerical 
results. 

Our model, initial conditions and simulation method are 
introduced in Section [2j The results of simulations are pre- 
sented and analyzed in Sections [3] and 2] Their comparison 
with a simple test model is discussed in Section [5] How 
can the results of the tilted stripping be estimated? This 
question is answered in Section [SI which deals with a mo- 
mentum delivery criterion. In the discussion of Section 
we compare our results with outcomes of previous numeri- 
cal simulations. 



2. Modeling 

2.1. Galaxy model 

O ur model galaxy is a late-type massive galaxy (LM-type 
in I Jachym et al.l I2007D : a three-component system with 
Plummcr bulge and dark halo, and with Toomre-Kuzmin- 
sech 2 (z/zo) disk. The values of the model parameters are 
M d = 8.6 10 10 M Q , a d = 4 kpc, z = 0.25 kpc, M b = 
1.3 10 10 M Q , a b = 0.4 kpc, and M h = 4.2 I0 11 M Q , a h = 20 
kpc. All three components are represented with particles 
(12 000 for the stellar disk, 12 000 for the halo, and 6000 
for the bulge). Their distribution is truncated at 16 kpc 
(disk) , 4 kpc (bulge) , and 40 kpc (halo) . The corresponding 
rotation curve is flat with a velocity of about 250 km s _1 . 
The ISM in the disk follows the same density profile as the 
stellar component with the same radial scale (ad — 4 kpc). 
Its mass forms 10% of the total disk mass. The ISM disk 
consists of 12 000 SPH particles which provides sufficient 
resolution for the large-scale effects of interest in this paper. 

2.2. Cluster model 

The galaxy cluster is modeled with dark matter (DM) 
and ICM gas. Their volume den s ities follow a /3-profil e 
(|Cavaliere fc Fusco-Femianol 119761: ISchindler et~atl Il999h : 
p = p (l + i? 2 /i? 2 )~ 3/3 / 2 . For the standard cluster model 
we use Picm = 1/2, pojcm = 6.5 10" 3 cm" 3 = 0.64 10~ 26 
g cm -3 and i? c ,iCM = 13.4 kpc. The fixed DM, which pro- 
vides the gravitational potential, has po,DM = 3.8 10~ 4 
M pc- 3 ,i? c ,DM = 320 kpc, and /3 D m = 1- The ICM is 
introduced with 120 000 SPH particles distributed within 
a 140 kpc truncation radius about the cluster center. 
Initially, the ICM particles are assigned velocity disper- 
sions calculated from the hydrostatic equilibrium equation 
of the cluster. The system then slightly relaxes to a sta- 
ble state. Particles that reach the truncation radius of 140 
kpc are relocated to the centrally symmetric position. See 
iJachvm et al.l (|2007[ ) for more details. 

2.3. Initial conditions 

The model galaxy initially starts with zero velocity from 
the cluster periphery (R = 1 Mpc) and freely falls along a 
strictly radial orbit towards the cluster center. In the stan- 
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Table 1. Set of basic simulation runs. We refer to Rlpl as 
to the standard run. Values of Eicm correspond to -Ricm c = 
140 kpc. 



run 


-Rc,ICM ; PO.ICM 


max 
rrajn 






(kpc,10~ 3 cm~ 3 ) 


cm _3 (km/s) 2 


(Mp,vc~ 2 ) 


R4pS 


53.6 , 52 


100 913 


114 




53.6 , 26 


46 752 


57 




OO.D , 0.0 


1U yyo 


1 i 

14 


m P o 


53.6 , 1.6 


2705 


3.6 


Rlp8 


13.4 , 52 


89 214 


41 


Rlpl 


13.4 , 26 


43 827 


20 


Rlpl 


13.4 , 6.5 


10 810 


5 


RlpO 


13.4 , 1.6 


2693 


1.3 


R0p8 


3.4 , 52 


86 549 


12 


R0p4 


3.4 , 26 


43 143 


5.9 


ROpl 


3.4 , 6.5 


10 769 


1.5 


ROpO 


3.4 , 1.6 


2689 


0.4 




T [Gyr] 

Fig. 1. Zero-order test vs. isolated evolution: Disk radii encom- 
passing 95 %, 90 %, 80 %, 10 % of the total ISM mass (from 
top down) as functions of time in the case of the galaxy evolving 
at rest in the center of the standard cluster (solid curves), and 
in an isolated galaxy (dashed). 



dard model the galaxy reaches the center at time T — 1.64 
Gyr with a velocity of about 1300 km s — 1 . The simulation 
is followed until T = 2 Gyr when the galaxy has already 
left the central part of the cluster filled with ICM particles. 
Thus, each simulation run corresponds to one orbit of the 
galaxy through the cluster center. 

To treat various cluster environments, we vary the val- 
ues of the i? c ,iCM and /?o,icm parameters, multiplying each 
standard value by factors of 8, 4, 1, and 0.25. Table [1] de- 
fines the basic set of simulation runs which are labeled with 
the corresponding factors (0 stands for the factor of 0.25). 
Table [1] gives for each run the combination of i? c ,iCM and 
Po.icm parameters, the peak value of the ram pressure along 
the orbit, and the column density of the ICM encountered 
on the orbit. Note that narrow ICM distributions or those 
with low values of density may represent ICM overdensi- 
ties or debris structures left over in the cluster from recent 
str ipping events . 

lAbadi et al.l (|1999D warned against too massive ICM 
particles that can punch holes in the gas disk and cause 
a large artificial drag on the ISM particles. In our stan- 
dard model (Rlpl) one ICM particle has a mass of about 
5 10 5 Mq, which is almost equal to the mass of an individ- 
ual ISM particle. Concerning the ICM component setup, 
its truncation radius of 140 kpc was chosen so that ram 
pressure in the standard cluster exceeds at that distance 
the gravitational restoring force in the disk outskirts of an 
approaching galaxy. 

We introduce the tilt angle of the disk with respect to 
the orbital direction of the galaxy as i = 0° for edge-on and 
i = 90° for face-on orientation, throughout the paper. 

2.4. Simulation method 

The simulations w ere carried ou t with the tree/SPH 
code GADGET JSpringel et alj 120011) adapted by 



Uachvm et ail (|2007h for calculations of ISM-IC M in- 
teractions. Smoothed particle hyd rodynamics (SPH, iLucvl 
Il977t iGingold k. Monaghanlll977l ) in its standard formula- 
tion has significant problems with contact discontinuities 
where the density jump is very large. Across such contact 
discontinuities, flui d instabilities can be suppressed (e.g. 
lAgertz et al.l 120071) . or othe r artifacts can occur when 
the sampling is poor (e.g. lOkamoto et alj 120031 ). The 



problem of ICM-ISM inte r action is basically a special 
case of this. Uachvm et all (|2007l ) modified GADGET to 
estimate smoothing lengths of either ICM or ISM particles 
separately from neighbors of a corresponding phase. 
This means that the search of the neighboring particles 
continues until a desired number of proper neighbors are 
found, while ignoring those of the second phase. It ensures 
that even a moderate number of ICM particles will fully 
cover during the interaction the ISM disk, as the ICM 
particles do not shrink when approaching t he dense ISM 
i n the disk. This modified code was used by I Jachym et alj 
(|2007l ) for simulations of face-on ram pressure stripping 
events. 

Although large ICM particles result in a low spatial res- 
olution and ISM particles that lack pressure gradients from 
the wind, our suggestion provides a simple way to avoid the 
above limitation of the standard SPH method, and it may 
be used for the quantitative purposes of this paper. An ideal 
solution could involve a new self-consistent way to estimate 
densities that allows a better representation of contact dis- 
continuities — e.g. based on a tessellation scheme (recently 
published by ISpringell2009D, implementation of the p arti- 
cle splitting mechanism (jKitsionas fc Whitworthll2007r i into 
the st andard SPH, or the new SPH scheme of Ka wata et alj 
(120091) . 

2.5. Zero-order simulation run 

In order to test the influence of the simulation technique 
used, and the effect of resolution on the results of our strip- 
ping simulations, we present a zero-order test: The galaxy is 
placed at rest at the center of the cluster and left to evolve 
surrounded by the ICM for about 300 Myr (the time com- 
parable to the duration of the standard ISM-ICM interac- 
tion discussed later) . The evolution of the galaxy's gaseous 
disk is displayed in Fig. [T] showing radii that encompass 
fixed fractions of the total ISM mass (solid curves) evolv- 
ing with time. The ISM disk (or at least its outer layers) is 
compressed by an abrupt introduction of the static external 
pressure of the ICM on a time scale of about 100 Myr. After 
this period the outer layers start to slightly oscillate while 
the disk slowly continues to settle towards an equilibrium. 
All the ISM particles during the simulation stay bound to 
the galaxy. The stellar disk component is unchanged dur- 
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Fig. 2. Pole-on views of the ISM disk crossing edge-on through the R4pl cluster. The disk rotates clockwise; the ram pressure 
operates from the left. Color saturation corresponds to increasing ISM density. 




T [Gyr] 

Fig. 3. Effect of the inclination: Evolution of the ISM mass 
bound to the galaxy (solid curves) or enclosed within an evalua- 
tion cylinder around the disk plane (dashed) in simulations with 
disk tilt decreasing from face-on to edge-on (from bottom up) of 
the standard Rlpl run. 



ing the simulation. For comparison, the evolution of the 
ISM disk of an isolated galaxy, not exposed to any external 
pressure, is plotted in Fig. [1] with dashed curves. 

The ISM component is treated isothermally in our sim- 
ulations. This simple approach is used to roughly mimic 
an efficient radiative cooling of the ISM gas. The disk con- 
finement in Fig. [T] may thus be understood in terms of an 
abrupt increase of the external pressure after the galaxy 
has been placed into the ICM - due to its isothermality the 
ISM cannot heat up (to compensate the increased external 
pressure) and is compressed instead. In Fig. [TJ an extreme 
effect is shown due to a maximum pressure corresponding 
to the cluster center. However, in the ram pressure strip- 
ping simulations that we deal with further in this paper, 
the galaxy arrives at the cluster center only gradually from 
its outskirts. This aspect will be treated in the following 
section. 

Concerning the influence of our numerical technique on 
the simulation outputs, no Rayleigh- Taylor instabilities or 
spurious effects due to the lack of resolution appeared in 
the zero-order test simulation. 



3. Process of stripping 

For non-face-on orientations of the galaxy with respect to 
the ICM wind, the galactic rotation starts to play a role. 
Its effect is most pronounced for near edge-on cases when it 
evokes an asymmetry in the gaseous disk. This can be seen 
in Fig. [2] where snapshots of the galaxy passing edge-on 
through the center of a wide RApl cluster are depicted. 
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Fig. 4. Mass flow of the ISM through the boundary of the eval- 
uation cylinder for five simulations of Fig. [3] 



(For illustrative purposes the RApl run is studied here 
rather than the standard one where edge-on stripping is 
less prominent.) Due to rotation, the disk side moving par- 
allel to the wind (upper side in Fig. ^ experiences ram 
pressure that is lower by a factor of (v — v rot ) 2 /(v + v TOt ) 2 
with respect to the opposite disk side; v is the orbital speed 
of the galaxy and u ro t is the rotational velocity of the disk. 
However, as e.g. the second panel of the Figure shows, the 
parallel-rotating side is more severely stripped than the op- 
posite one and it is thus not the highest ram pressure that 
yields the strongest strippi ng. This paradox was noted by 
iRoediger fc Briiggenl (|2006l ) as well. The explanation is in 
terms of time delay: on the parallel-rotating side, the gas 
is accelerated towards more rotational velocity and exits 
the galaxy almost at once. At the opposite side, the gas is 
being braked from v — v mt towards v and with lower or no 
centrifugal force it has time to fall towards the center of 
galaxy. Then, it still has to rotate half a turn around the 
galaxy, and it has time to be re-accreted (or at least it stays 
around the galaxy for a longer time). Snapshots in Fig. [2] 
further show how, due to rotation, the tail of the stripped 
material winds up around the disk. 

To measure the actual level of stripping we set up a 
cylindrical zone about the galactic disk plane (r = 16 kpc 
and \z\ < 3 kpc), and follow the enclosed ISM mass and 
the flow of the material through the zone boundary. This 
is shown in Figs. [3] and 0] for five simulations of the stan- 
dard run with different tilt angles, together with the ISM 
mass fraction bound to the galaxy (solid curves in Fig. [3j . It 
confirms the expected trend: with the tilt angle growing to- 
wards face-on, more material is stripped. The maximum de- 
parture of the curves of the bound and in-the-zone-enclosed 
ISM mass in Fig. [3] indicates the amount of material that 



Jachym et al.: Ram pressure stripping of tilted galaxies 



5 




1.5 1.6 1.7 1.8 1.9 2 



T [Gyr] 

Fig. 5. ISM re-accretion: Distances behind the galaxy (in the 
direction of the wind) within which a certain fraction of actually 
bound ISM particles are enclosed (from 95 % to 80 %, curves 
from top down). Face-on (black) and edge-on simulations of Rlpl 
run. (In the edge-on case the distances are taken in absolute 
value as initially one half of the particles occur ahead of the 
galactic mid-plane.) 



is re-accreted after the ISM-ICM interaction. The bound 
mass in Fig. [3] is given by the sum of the ISM particles for 
which E^in + E pot < 0. Note that the system of the galaxy 
is not closed - it experiences a changing effect of the ex- 
ternal potential of the cluster and an external pressure. A 
slight growth of the solid curves in Fig. [3] may be found at 
later simulation times. 

The "stripping rate" M TSM , i.e. the flow of the ISM 
through the boundary of the evaluation zone, exceeds al- 
most 400 M©yr _1 for a face-on galaxy, the peak value 
decreases towards the edge-on orientation to about 50 
M Q yr _1 (see Fig. 0}. This Figure shows not only the re- 
moval rate (Mism < 0) but also a reversed re-accretion rate 
of the shifted material (Mism > 0) - in the face-on case it 
reaches about 50 M Q yr _1 after the passage of the galaxy 
through the central cluster part. (Note that the amount of 
re-accretion is dependent on the extent of the evaluation 
zone.) 

The re-accretion process may be better followed in Fig.0 
which traces the actual distances of the bound ISM behind 
the galaxy (in the direction of the wind) for face-on and 
edge-on cases. In the face-on orientation, the ISM is by an 
increasing ram pressure first gradually pushed out of the 
disk plane, and then rapidly "thrown" to large distances, 
shortly after the cluster center passage. This material re- 
turns to the disk at later simulation times. In the edge-on 
orientation, the disk is first compressed by the rising ram 
pressure, and then some material is shifted to larger dis- 
tances. The peak in Fig. [5] is, however, much smaller than 
in the face-on case. The reason why in near edge-on orien- 
tations the pushed material stays close to the disk is in the 
conservation of the angular momentum of the ISM parti- 
cles rotating in the disk. When pushed by the ram pressure, 
the particles wind up around the disk. We checked this be- 
havior in a more efficient run R4pl where almost the same 
amount of ISM is stripped as in the standard face-on run. 

Apart from the total stripped amount, no radical differ- 
ences in the duration of the stripping process of variously 
tilted galaxies can be seen from our simulations. Several de- 
tails, like a slight time shift of the maximum of the stripping 
rate, or the course of re-accretion, look different in edge-on 




0.1 0.2 0.3 



T [Gyr] 

Fig. 6. Ram pressure stripping vs. zero-order test: Actual esti- 
mate of the stripping radius of the face-on galaxy in the standard 
Rlpl cluster (solid) as a function of time, vs. radius of the disk 
confined by static ICM pressure (dashed, the 95 % radius from 

Fig.ru. 

and face-on orientations, however these effects seem to be 
of secondary nature. 

3.1. Static vs. dynamic ICM pressure 

How does the effect of the disk radius reduction due to 
ram pressure compete with the disk confinement due to the 
static ICM pressure discussed in Section [2T5T ? In Fig. [6] the 
rate of the disk radius reduction caused by the ram pressure 
stripping in the standard face-on run is compared with the 
upper curve of Fig. [TJ It shows that the effect of the ram 
pressure is more important, even though the zero-order test 
took place in the center of the standard cluster where the 
disk experienced the maximum external pressure possible. 
In the ram pressure simulation the in-flying galaxy enters 
the ICM distribution at the truncation radius where the 
ICM density (and its pressure) is only small, and plunges 
deeper towards the cluster center, which is crossed very 
fast. Thus the actual role of the static ICM pressure will be 
even less than in Fig. [SJ 

Although the disk confinement due to the static ICM 
pressure operates against ram pressure - it does not strip 
the gas but makes the disk more resistant to the ram pres- 
sure - the tendency of the increasing effect of the static 
ICM pressure in "larger" clusters goes in the same di- 
rection as the effec t of ra m pressure stripping. Recently, 
iTonnesen fc Brvar] (|2009l ) have taken into account the 
multi-phase nature of the gas disk and found that the strip- 
ping radius is considerably smaller for galaxies with radia- 
tive cooling, although the amount of stripping is similar 
to that in galaxies without cooling. This goes in the same 
direction as the above effect observed in our simulations. 

We conclude that the disk confinement due to static ex- 
ternal ICM pressure is much less important than the strip- 
ping due to ram pressure and we neglect this effect in ana- 
lyzing the simulations in the following sections. 

4. Results 

In this section the results of the basic set of simulation 
runs introduced in Table [T] are described. For each run 
listed we have performed five simulations with increasing 
disk inclination: i = 0° (edge-on), 20°, 45°, 70° and 90° 
(face-on). The overall Table provides the results: mass 
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Table 2. List of performed simulations - results. From left to right: fraction of the stripped ISM mass for simulations with 
decreasing inclination, stripped mass fraction estimate of GG72 criterion, corresponding stripping radii and the estimate of GG72. 



run 


M 9 ?° 

-"-'strip 

(%) 


M 7 t °° 

-"-'strip 

(%) 


(%) 


(%) 


■""strip 
(%) 


M^ 72 
(%) 


r 90° 
strip 

(kpc) 


r 70° 
strip 

(kpc) 


r 45 ° 

strip 

(kpc) 


„20° 
strip 

(kpc) 


r° ■ 

' strip 

(kpc) 


r GG72 
strip 

(kpc) 


R4pS 


93 


91 


89 


86 


70 


93 


1.5 


0.9 


1.4 
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fraction of the stripped material (Mgtrip) at the end of the 
simulation for five disk inclinations, stripped mass fraction 
estimated with the GG72 formula, stripping radii r stri p for 
five inclinations, and the stripping radius corresponding to 
GG72. 

4.1. Stripped amount 

The fraction of the ISM particles that is no longer bound to 
the galaxy at the final simulation time yields the stripped 
mass fractions, whose dependence on inclination is shown in 
Fig. [7] The general tendency of the plot is evident - with the 
inclination decreasing from face-on to edge-on, the amount 
of stripping declines. Furthermore: 

— Previous simulatio ns of iQuilis et alj (|2000j) . 

iRoediger fc Bruggenl (|2006l ) and others have shown 
that there is almost no difference between i = 70° and 
face-on stripping (i — 90°). This is true in Fig. [7] for all 
runs. 

— For large pressure peaks (runs i?4p4 and RAp8), the 
stripping amount is almost independent of the inclina- 
tion (for angles higher than about 20°). 

— The dependence of the stripping amount on inclination 
is more pronounced for smaller ram pressure peaks. 

— Runs with the same value of i? c ,iCM ■ Pojcm quantity 
show close profiles of the M s t r i p (i) curves - e.g. (i?4pl 
and RlfA), or {R4p0, Rlpl, and R0p4). 

The behavior with varying inclination is more clearly 
displayed in Fig. [5] which shows a new parameter r\ that we 
name the stripping efficiency with respect to inclination. It 
is the ratio of the stripped mass fraction relative to the 
face-on case: 

,.s Afjjtrip.i 1 — Mfl n i/MixA ,^ 

Tj[t) = — — = ■ L - L -r . (1) 

•M s tri P) go° 1 — Mfi ni goo/Mi n i 

Thus, n characterizes the relative strength of a given ram 
pressure profile to strip the ISM from an inclined galaxy 
with respect to the face-on case. In Fig. HI one observes 
that the stripping efficiency always declines with decreas- 
ing inclination angle. Moreover, both wider and higher ram 
pressure peaks yield higher efficiencies. For smaller pressure 
peaks, rj falls fastest at medium inclinations (between 20° 
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Fig. 7. Stripped mass fraction as a function of disk inclination 
for runs from Table [2] Each run is labeled with the peak value 
of the ram pressure. The inclination is 0° for edge-on and 90° 
for face-on orientation. 




inclination 



Fig. 8. Stripping efficiency with respect to inclination n — 
M s trip,i/M a trip,9o° as a function of the inclination (runs i?0pl 
and .ROpO are not plotted). 

and 70°), while for large peaks at inclinations smaller than 
20°. 

Stripping of face-on galaxies depends on the amount 
of the ICM encoun tered along the trajectory (Fig. 15 in 
iJachvm et al.l I2007T ): galaxies crossing clusters with such 
combinations of i? c ,icM and Po,icm parameters that Sicm ~ 
const were stripped to the same level. For inclined cases 
this hypothesis is tested in Figs. r9land[TUl which show the 
stripped mass fraction and the stripping efficiency with re- 
spect to inclination as a function of the column density of 
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Fig. 9. Stripped mass fraction as a function of the column 
density of the encountered ICM. The symbols denote the re- 
sults of the SPH simulations for various inclinations. The curves 
correspond to the fitting curve in Eq. [2] 



the encountered ICM, respectively. Different colors refer to 
simulations with different inclination angles. The following 
trends are apparent 

— with increasing amount of encountered ICM the 
stripped mass fraction and the efficiency increase; 

— there is a tendency for the results from each inclina- 
tion to follow a certain close relationship, despite some 
systematic deviations due to the different ICM density 
distributions; 

— for high SicMj these relations saturate towards com- 
plete stripping; 

— for lower ICM column densities, edge-on stripping is re- 
duced with respect to face-on by a constant factor. 

All these results indicate that the column density of 
the ICM encountered during the ISM-ICM interaction is 
the key parameter determining the stripping outcome. It is 
much more important than the maximum value of the ram 
pressure experienced along the orbit (GG72 criterion). 

The numerical results in Fig. [5] can be fitted with the 
following formula 



AL 



strip 



Mi, 



I-exp [-0.01E IC m(1 



6sin 15 i)] . 



(2) 



The exponent 1.5 provides a rather satisfactory match to 
the results with intermediate inclinations, as can be seen in 
both Figs. M and [TO] 

In Table [2J the stripped mass fraction estimated with 
the GG72 criterion is added. For the largest ram pressure 
peak (i?4p8), the values match well the face-on results of the 
SPH simulations. For smaller peaks, however, GG72 overes- 
timates the real stripped fraction. Of course, no dependence 
on the tilt angle is included in the GG72 criterion, and its 
multiplication by a simple cosi term does not resolve the 
problem, as it yields zero stripping in edge-on cases. 

4.2. Stripping radii 

Due to an asymmetric shape of the gas disk that develops, 
especially after an inclined stripping event, the determina- 
tion ofiteffnalrax^ radius, is not trivial 
(see iRoediger fc Bruggenl l2b06'). In close to edge-on cases, 
all the stripped and shifted material stays in the plane of 
the disk and makes the determination difficult. 
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Fig. 10. Stripping efficiency with respect to inclination as a 
function of the column density of the encountered ICM. The 
curves correspond to the fitting curve in Eq. [2] 




Fig. 11. Stripping radii r atr i p versus column density of the 
encountered ICM. The symbols are the same as in Fig. [9] The 
points are fitted with the fitting function in Eq. [3] 



In Table[2j the given values of stripping radii correspond 
to radii within which 90% of the bound ISM is enclosed in 
the final simulation time. This may slightly underestimate 
the radius for the face-on orientation, while overestimating 
it for the edge-on case. The values do not show a strong de- 
pendence for large and small peak ram pressures; however, 
for medium pressure the trend of increasing stripping radii 
towards edge-on orientations is visible. 

As shown in Fig. Ill) the stripping radii for all simula- 
tions can be quite reasonably fitted with 



r str 



I2kpc / [1 + Eicm(0.1 + 0.15 sin i)}. 



(3) 



The stripping radii estimated in the last column of 
Table [2] using the GG72 criterion show quite a good agree- 
ment for the largest run. However, there is a tendency to 
overestimate the SPH results. For smaller peaks, GG72 pre- 
dicts less stripping corresponding only to the peak ram 
pressure. 

5. A simple test particle model 

In order to study the influence of the parameters more effi- 
ciently, and to separate various effects - the direct kinemat- 
ics on the ISM induced by the ram-pressure, and the related 
hydrodynamical consequences - we construct a very simpli- 
fied model for the removal of gas from a galaxy. The gas 
disk is represented by 10 000 point masses distributed as in 
an exponential disk and held in the potential well composed 
of the bulge, stellar disk, and dark halo, with parameters 
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identical to those used in the SPH computations. Initially, 
all these particles are in circular orbits about the center of 
the galaxy. Then, they are subjected to a Gaussian-shaped 
force pulse of specified tilt angle with respect to the galactic 
plane, duration, and maximum amplitude. The simulation 
starts at a time that is three times the pulse duration before 
the maximum. In the computation of the acceleration, the 
inertial mass of each particle is taken to be proportional 
to the gas surface density of the exponential disk, in order 
to account for the correct computation of the ram pressure 
in that part of the disk. For each particle, the equation of 
motion is integrated, and when the total energy becomes 
positive and the particle becomes unbound, it is marked 
and counted as having been stripped. If all particles were 
to experience this force pulse, the fraction of stripped parti- 
cles is independent of the tilt angle, because all reach escape 
velocity within time-spans that depend little on their initial 
position in the galaxy and the force tilt angle. 

In order to simulate the thickness of the disk, we in- 
troduce a shielding or shadowing of the particles down- 
stream of those particles facing the f orce. In this, w e follow 
the modeling by sticky particles by IVollmer et al.1 ((2001) . 
Particles that are shielded by at least one upstream particle 
are assumed to experience only the galaxy's internal forces. 
Experiments show that with the choice of a suitable value 
for the shadowing radius (typically 200 pc), we are able 
to reproduce well the amount of stripping for the standard 
model of the SPH computations, as well as its dependence 
on the tilt angle. 

As is seen in Fig. [12J the stripped mass fractions from 
these simplified models behave in a very similar way as the 
SPH results. While at high values of the ICM column den- 
sity the stripping is less severe and the relation is slightly 
flatter, the influence on the tilt angle is nearly identical 
- one obtains the same amount of stripping in the edge- 
on configuration as in the face-on case for an ICM column 
density that is larger by about a factor of 6. This indicates 
that this modeling appr oach, and hen c e also the sticky- 
particle computations by IVollmer et al.l (|200lh , appears to 
be quite reliable in capturing the main features of ram- 
pressure stripping. However, this simplified approach has 
notable differences to the full hydrodynamical approach, in 
features that we consider below. 

The stripping efficiencies with respect to inclination 
(Fig [r3]) show that the ratios of face-on and edge-on re- 
sults are as large as from the SPH models, but that there 
are some limitations. First of all, at intermediate angles the 
simplified modeling gives results that are closer to the face- 
on case than we find in SPH. Furthermore, at high ICM 
column densities, one does not find such a steep increase 
of stripping as experienced in SPH. We searched quite ex- 
tensively in parameter space, but could not find situations 
where the behavior matched the SPH results more closely. 
Evidently, this is a consequence of the limitations of the 
simplified treatment. We take this as an indication that 
this kind of diagram is sensitive to the different treatment 
of the (hydro-)dynamics of the problem. 

The simplified modeling does allow us to separate the 
different physical phenomena in the ram-pressure process. 
In particular it helps to identify the hydrodynamical shield- 
ing by the matter lying upstream as the principal cause for 
the dependence of stripping on tilt angle, seen in our SPH 
results. The deviations seen with respect to the SPH results 
suggest that the details of the behavior changes with tilt 
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Fig. 12. The dependence of the stripped mass fraction on the 
column density of the swept-up ICM, for different tilt angles. 
The dots are results from a grid of simplified point-mass models, 
covering a range from 1/16 to 16 times in both duration and 
strength of the ram pressure force pulse, centered about the 
value for our standard model. The curves are the analytical fit 
to the SPH results. The symbols and colors are identical to those 
of Fig. [9] Note that the range of parameters in both axes is much 
wider. 
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Fig. 13. Similar to Fig. 1121 but for the dependence of the strip- 
ping efficiencies with respect to inclination as a function of the 
column density of the encountered ICM. 



angle comes from the finer details of the hydrodynamics 
involved in the problem. 



6. Comparison with the analytical criterion 

From the consideration of the momentum provided by 
the in-stream i ng IC M to the gas in the galactic disk 
I Jachym et al.l (|2007t ) derived the condition under which a 
volume element of disk gas (with column density Eism for 
the face-on situation) would escape from a galaxy, 



(v) 



Picm^ICM 



ISM 



Rafter > ^csc = \/-2$(r), 



(4) 



i.e. if the speed v a fter exceeds the local escape speed. The 
stripping radius would then be defined as the radius where 
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Rafter = ^osc and the stripped mass fraction as the mass 
exterior to that radius: 



-^strip — 



2n^ lSM {r)rdr / M{ 

' out 1 ) 



(5) 



where r out is the truncation radius of the disk. For an 
exponential disk EismM = So exp(— r/a,d) with A/ism = 
27rSo a d w e have 

(o d + ra„t) e" r °-/ a « - (a d + r strip ) e -»W«* 
JWstrip — -, ■ \ — -T- ■ Id J 

{ad + r out ) e r »»'/"" - a d 

On the other hand, the escape velocity as a function of 
galactocentric radius is the sum of the contributions from 
all components, i.e. bulge, disk, and dark halo. At the 
stripping radius the escape velocity equals the speed w a fter, 
which is a function of the ICM column density which we 
shall now take as the independent variable. Hence we can 
compute the stripping radius as a function of ICM column 
density, and with that we get the stripped mass fraction. 
This is easily evaluated numerically by computing for each 
galactocentric radius the mass fraction outside of that ra- 
dius, and on the other hand the quantity 



Sicm = 



-aSM 



(r) 



-"ISM 



(7) 



( V )piCM ( V )piCM 

Figure [HI depicts the relation between M str i p and Eicm for 
our standard model. The above equation is able to match 
the SPH results very closely, if we adopt A = 0.5. In this 
diagram, our fit formula to the SPH results (Eq. ^ gives a 
function very similar in character to the one derived from 
our analytical criterion. The choice of the outer disk radius 
influences the shape of the curve in Fig. [JJ] we find that a 
value of 20 kpc gives a better match than the value of 16 
kpc assumed for the SPH models. This detail is not incor- 
porated in the fit formula, whose results therefore match 
the SPH results at low stripping less well. 

Although the SPH results from the edge-on case do ex- 
hibit a larger scatter than the face-on data, it is quite re- 
markable that they can still be matched by the analytical 
prediction for the face-on case, if one simply increases the 
necessary ICM column density. This means that the edge- 
on models can well be considered as having a lower stripping 
efficiency, with a reduction factor of about 5. This factor 
may be regarded as depending only on the tilt angle i, as 



A(i) = 0.12 + 0.38 sin i. 



(8) 



A sine-type dependence could be understandable as a ge- 
ometrical projection effect, in that the edge-on galaxy 
presents a smaller cross section to the ICM. If it were such a 
pure effect of the effective cross sections, one should expect 
for our standard model a ratio of the areas of 4/0.25 = 16, 
from the disk radial scale being 4 kpc and its height being 
0.25 kpc. For the moment, we prefer to take our analyti- 
cal formula as a convenient expression which encapsulates 
essential features of the problem rather than as giving a 
complete description. 

7. Discussion 

We compare our results to the findings of previous numer- 
ical studies which treated inclined ram pressure stripping. 
These investigations differ in their numerical treatment of 
hydrodynamics, the model for the galaxy, and the way the 
ICM gas flows around the galaxy. 
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Fig. 14. The relation between stripped mass fraction and the 
column density of the swept-up ICM gas, computed for our stan- 
dard model from our analytical criterion. The dots are results 
from the SPH models for the face-on configuration, compared to 
the predictions, adopting a value of A = 0.5 {solid curve). The 
open circles are results in edge-on orientation, and the dashed 
curves are the predictions with an adopted value A — 0.12. The 
dotted curves come from the fitting formula Eq. [2] 
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Fig. 15. Str ipped mass fraction as a func tion of the disk incli- 
nation , from IVollmer et "ahl (|2001 | ) - bl ack, [Roediger & Briiggcn 
(120061) - red IRoediger fc Brugger] (|2007h - green, and 
ISchulz fe Struckl <|2001l ) - blue. Labels giv e values of the peak 
ram p ressure, in units of cm _3 (km/s) 2 . For Roe diger fc Briiggenl 
(2007) the actual values of the inclination are given for the in- 
stant of pericenter passage. 



7.1. General trends 

The fra ction of the ISM mass that is stripped in the simula - 
tions of [Vollmer et al.l (120011) . [Roediger fc Briiggenl (120061) . 
IRoediger fc Briiggenl (|2007t >. and lSchulz fc Struckl (|200ll ) is 
plotted in Fig. [15] as a function of inclination angle. The 
curves join runs with different disk inclinations, and are la- 
beled by the maximum ram pressure. It must be stressed 
that for different curves not only does the maximum value 
of the ram pressure vary, but also its duration. These gen- 
eral features are apparent in Fig. 1151 



for inclinations larger than about 50° there is very little 

dependence of the stripped fraction; 

towards the edge-on orientation, on the contrary, the 

amount of stripping drops strongly; 

the curves flatten for higher peak ram pressures. 
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Fig. 16. Stripping efficiency with respect to inclination as a 
function of the angle for the same results as in Fig. 1151 



Our results show the same features (cf. Fig. [7]). A 
clo ser comparison show s that the near edge-on models 
by IVollmer et al.1 (|2001h are substantially more strongly 
stripped than ours for the same peak ram pressure. 
However, since their Fig. 7 suggests that the radial scale- 
length of the gas disk is longer than in our model, one 
should expect that due to the higher ISM surface density 
they should be stripped less s everely. This suggest s that the 
shielding mechanism used bv lVollmer et al.l (|200l| ) is overly 
effective in near edge-on orientations and thus overesti- 
mates the stripping, at le ast in comparison to t he fu ll hy- 
drodynamic treatment of iRoediger fc Briiggcn (2006) and 
our SPH approach. 

For the same values of the stripped mas s fraction, 
the simulations by IRoediger fc Briiggenl (|2006D require a 
substantially lower value of the peak ram pres sure than 
our models or thos e by Vollmer et al.l (|2001[ ), because 
IRoediger fc Briiggenl (|2006 ) consider the flow of the ICM 
to be constant in time, hence much longer than in the other 
models. 

The stripping efficiency with respect to inclination (rj) 
measures the relative stripping potential of individual incli- 
nations with respect to the face-on case. It might thus serve 
as a good instrument for comparing results from different 
numerical codes. From Fig. [11] it is obvious that the results 
of different authors give very similar dependences on the 
inclination angle. In spite of the basic difference in the ab- 
solut e values for the maximum r am p ressure, the results o f 
both IRoediger fc Briiggenl (|2006t) and I Vollmer et al.l (|200lD 
show that higher ram pressures result in a lower depen- 
dence of the stripping on inclination angle, i.e. making the 
inclined stripping as effective as the face-on case. Thus, the 
features of the inclination dependence, which our present 
study shows in a systematic way, has already been seen in 
earlier studies. 



7.2. Vollmer's fitting formula 

Voll mer et al.l (|2001[) derived from their simulation results a 
formula which relates the fraction of the final and initial to- 
tal gas mass to the peak ram pressure and inclination angle 
(their Eq. 21): M ini /M Rn = 0.25p™Z sin 2 [0.9( i + 10°)] + 
0.84, where 1S normalized by po — 100 cm~ 3 (km/s) 2 . 
Using this fitting formula we reproduce Vollmer's Fig. 17 
in Fig. [T7] as a function of angle we plot the stripping effi- 
ciency with respect to inclination r\. It shows that although 
the fitting function does not match the simulations per- 
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Fig. 17. Curves of stripping efficiency with respect to inclina- 
tion rj a s a function of tilt ang le as calculated from the fitting for- 
mula of lVollmer et al.1 (|200ll ) for various values of the peak ram 
pressure: 100, 300, 1000, 2000, 5000 and 10 000 cm" 3 (km/s) 2 
(from left to right). Results of Vollmer's simulations (dots) cor- 
respond to the solid curves. Results of our i?4p0 and R4pl runs 
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function of the column density of the encountered ICM for simu- 
lation s of lRoediger fc Briiggenl l|2006f ) (green) and IVollmer et all 
(|2001h (black). 



fectly (and even gives for low ram pressures negative values 
of rj) the general trend corresponds well to the numerical 
results, and to our findings as well: that for high ram pres- 
sures the efficiency of stripping is high and less dependent 
on the tilt angle than for low ram pressure values. 

7.3. Dependence on ICM column density 

To look for the dependence of the stripping results on 
Sicm in the previous simulations, we estimate these values 
from t he parameters of their runs. In Roedi ger fc Briiggenl 
(200(f), a constant flow of ICM with three values of den- 
sity (0.01, 0.1, and 1 x 10 _26 gcm~ 3 ) provokes a wind 
at 800 kms -1 lasting for 1 Myr. This yields values of 
the ICM column density of 0.4, 4, and 36 M pc~ 2 for 
their weak, medium , and strong pressure models, respec- 
tively. IVollmer et al.l (|200l[ ) approximate ram pressure pro- 
files with Lorentz functions, whose duration is determined 
by the condition that 500 Myr before the peak, a constant 
pressure of 50 cm™ 3 (km/s) 2 is reached. For their four peak 
ram pressures of 1000, 2000, 5000, and 10 000 c m - 3 (km/s) 2 
we derive widths £hw = 115, 80, 50, and 35 kpc, respec- 
tively, which yield values of Sicm: 4.4, 6.2, 9.4, and 13.5 
M Q pc- 2 . 
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Figure [TH] compares the t?(Sicm) d epen- 
dence for the res ults of IVollmer et al.l ( 20011 ) and 

iRoediger fe Briiggenl (120061). We do not include the 
results of iRoediger fc Bruggenl (|2007t ). as the tilt angle 
changes along the g alaxy's orbit. The results of both 
IVollmer et alT (|200lli and IRoediger fc Bruggenl (|2006f) 
show the same trend as that revealed in our simulations: 
with increasing Sicm the inclination efficiency grows. 
It is worth noti ng that the increa se is much steeper 
in the case of IVollmer et al.l (1200 lh . in comparison to 
IRoediger fc Briiggenl (|2006l ) and our own results. Since the 
parameters of the model galaxies are rather similar and 
the stripping efficiency r\ measures the relative potential of 
inclined stripping with respect to the face-on case, it could 
well be that these differences are due to the different ways 
the hydrodynamics is treated. 



7.4. Limitations of strictly radial orbits 

To separate the effects of the disk tilt angle on the stripping 
results we have focused in this paper on galaxies on strictly 
radial orbits. Then the tilt angle stays constant along the 
orbit. In more realistic orbits however, the disk inclination 
usually changes with time. Thus, almost all galaxies are 
likely to experience once along their orbit an efficient (close- 
to) face-on ram pressure phase. Its duration and strength is 
then determined by the parameters of the orbit. As the re- 
sults of our simulations have shown, the tilt-variance might 
be important for situations with compact ICM distributions 
where the face-on stripping really is more efficient than the 
tilted ones. In extended clusters, however, the dependence 
on the disk inclination vanishes as the stripping tends to 
saturate. 

Concerning the stripping potential of individual models, 
note that the strictly radial orbits used in our calculations 
may be considered as being equivalent to slightly elliptical 
ones: The ram pressure profile along a strictly radial orbit 
through a given ICM distribution may well be the same as 
the ram pressure profile along an elliptical orbit, however 
through a modified ICM distribution that is more centrally 
peaked (higher) and less extended (more compact). In other 
words, on an elliptical orbit the galaxy passes through the 
slope of the ICM density peak. 



8. Conclusion 

We have studied in detail the influence of disk tilt angle 
on the ram pressure stripping efficiency in galaxy clusters. 
A grid of N-body/SPH numerical simulations, varying the 
density of the ICM, the peak ram-pressure, its duration, 
and the orientation of the galaxy was performed. We sum- 
marize below the main results of the simulations: 

— the stripping amount is lower for galaxies moving edge- 
on with respect to the ICM wind, in agreement with 
previous simulations; 

— apart from the total stripped amount, there is no radical 
difference in the duration of the stripping process of 
variously tilted galaxies; 

— due to conservation of angular momentum, ram pressure 
shifts the ISM in near-edge-on orientations to smaller 
distances behind the galaxy than in the face-on case; 



— the dependence of the stripping amount on the disk tilt 
angle vanishes for large clusters, where the amount of 
the encountered ICM (Sicm) is large; 

— on the other hand, for galaxies crossing small ICM over- 
densities or debris structures left over in the cluster from 
recent stripping events, the dependence on the tilt angle 
is more pronounced; 

— we emphasize the role of the amount of the encountered 
ICM to the ram pressure stripping results; 

— the stripping efficiency with respect to inclination rj (i.e. 
the ratio of the stripped mass fraction relative to the 
face-on case) depends strongly on the disk tilt angle for 
small ram pressure peaks but again vanishes for large 
values of Eicm when stripping saturates; 

— for non-face-on orientations, the sense of rotation plays 
a role: contrary to expectations, it is not the side with 
highest ram pressure which is more stripped; 

— we propose a fitting formula (Eq. [2]), which should now 
be tested against observations, deriving the ICM density 
from X-ray imaging of clusters, and the gas stripping 
from HI observations. 

By comparison with a simple analytical stripping con- 
dition (| Jachym et al.ll2007t ). the total momentum imparted 
by the encountered ICM on the galaxy gas appears to be 
the main determinant of the stripping resul t. Contrary to 
the G G72 criterion or the fitting formula of IVollmer et al.l 
(1200 lh who assign the key role to the peak value of the 
ram pressure, we stress the importance of the dependence 
on the amount of the encountered ICM, together with the 
disk inclination angle: 

Mstrip/Mnit ~ SicM(a + b sin Q i). 

The values of a, 6, and a depend on the ISM distribution in 
the unperturbed disk. We use such a dependence to fit the 
outcome of our numerical simulations - the stripped mass 
fraction and the stripping radius. The analytical stripping 
condition based on the momentum balance yields similar 
formulae. 

Using simple test particle models taking schematically 
into account the "shadowing" of the gas disk by the 
material occurring upstream, we can roughly obtain the 
main features of the stripping. The imparted momentum 
("^picm^icm proves to be the key parameter of the ICM- 
ISM interaction. Such particle models however fail to reflect 
the true dependence of the stripping efficiency on ICM den- 
sity and inclination revealed by the full hydrodynamics. 
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